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I. INTRODUCTION 

When a homogenous system departs suddenly from equilibrium, it will spontaneously 
segregate into two different states, thermodynamically more stable. 

This process can either initiate via a nucleation process, where an energy barrier has to 
be crossed, or via a spinodal decomposition when the system is led into a linearly unstable 
configuration. In this latter case, the leading instability selects a modulation of the order 
parameter at a well defined length scale. This instability will grows and, due to nonlinearity, 
saturates. The resulting micro-segregated pattern is composed of well defined interfaces 
delimiting monophasic domains containing one of the two stable phases. These interfaces 
will then interact with each other and coalesce, during a much slower, self-inhibiting process, 
where the number of domains diminishes whereas their typical size increases, ending with the 
formation of either a micro-segregated complex pattern or a single interface separating 
two semi infinite domains, one for each new stable phase (macro-segregation). 

Hillert^, Cahn and Hilliard[3?] have proposed a model equation describing the segregation 
for a binary mixture. This model, known as the Cahn-Hilliard equation (C-H later on), 
belongs to the Model B class in Hohenberg and Halperin's classification [4]. It is a standard 
model for phase transition with conserved quantities and has applications tophase transition 
in alloys [2|, binary mixtures [5|], vapor condensation 6[] and liquid crystals [TJ, segregation of 
granular mixtures in a rotating drum[s], or formation of sand ripples jsl, [lo| . 

In this article, by spinodal decomposition, we refer to the first stage of the dynamics only, 
while coarsening or Ostwald ripening will denote the second stage. Although this coarsening 
dynamics is in fact already present, its influence can be neglected during the first stage of 
the process. 

An important activity has been devoted to the description of the dynamics of phase 



transition, using experiments, numerical simulations [11] or scaling methods [12J. The late 
stage of the dynamics, where the Ostwald ripening dominates, exhibits "dynamical scaling" : 
the dynamics presents a self-similar evolution where time enters only through a length scale 
L(t), associated with a typical length of the domains. This scaling argument gives the law 
Lit) ~ t 1//3 for spatial dimensions greater than one and a logarithmic behavior in D=l. 



This last stage, as observed in two-dimensional demixion of copolymers [l^] and as sug- 



gested initially by Langer 



14j . can be described as a process of synchronous fusion and 



evaporation of domains, spatially alternated. It corresponds to the evolution that breaks 
the least the symmetries of the pattern. This observation motivated our work and the aim 
of this article is to make use of a one dimensional ansatz to describe quantitatively the ideal 
coalescence process. This ansatz is in the form of a one parameter family of symmetric pro- 
files which interpolates between two stationary states composed of homogeneous domains 
of length A and 2A. It allows in ID to describe a self similar sequence of coalescence and 
evaporation of domains (or of their dual counterparts, the interfaces) , starting from the pe- 
riodic micro-segregated state which ends the spinodal decomposition dynamics, and leading 
continuously to either a single interface, or a stable ID modulated phase. 

The paper is organized as followed: part 2 will focus on the first part of the dynamics, 
i.e. spinodal decomposition. We will reproduce briefly the original derivation by Cahn 
and Hilliard, restricting ourselves to the one dimensional case, mainly to fix the notations 
(part 2.1). In part 2.2, a family of symmetric solutions of the Ginzburg-Landau equation 
will be used to study the non linear part of dynamics and in part 2.3 we will identify all 
the intermediate symmetric stationary states of the (C-H) dynamics. Then in part 3, we 
will turn to Ostwald ripening : in part 3.1, a non-symmetric family of solutions of the 
(G-L) equation is used to construct a continuous interpolation between two consecutive 
symmetric stationary states. After a study of the energy landscape associated with this 
ansatz (part 3.2), we will compute in part 3.3 the characteristic time associated with one 
step of coalescence : we will recover the ID logarithmic law. Before concluding, as an 
application in part 4, our ansatz will be used, in the case of micro-segregation, to compute 
the period of the final thermodynamically stable modulated pattern. 
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II. THE CAHN-HILLIARD MODEL 



A. Linear stability Analysis 

The Cahn-Hilliard (or Conservative Time Dependant Ginzburg Landau) equation is a 
modified diffusion equation for the scalar order parameter which reads in its dimensionless 
form: 

f (r,i) = V 2 ^^ = V 2 (|« + 2^-V 2 *). (1) 

The real order parameter can correspond to the dimensionless magnetization in Ising fer- 
romagnet, to the fluctuation of density of a fluid around its mean value during a phase 
separation or to the local concentration of one of the components of a binary solution, e is 
the dimensionless control parameter of the system ; it is often identified to the reduced tem- 
perature (e = T ~ Te where T c is the critical temperature of the first order phase transition). 



This equation, first derived by Cahn and Hilliard, has also been retrieved by Langer(14| 
from microscopic considerations. A conservative noise can be added to account for thermal 
fluctuations [15(, but in this article, we will only consider the noiseless (C-H) equation. 

It admits homogeneous stationary solutions which are extrema of the symmetric Landau 
potential V(\&) = l^ 2 + |^ 4 - For positive e, there is only one homogenous solution \& = 
which is linearly stable. A pitchfork bifurcation can be experienced when quenching the 
system from a positive reduced temperature e to a negative one : the — solution 
becomes unstable ; two other symmetric stable solutions appear ^ = ±^p. 

Spinodal decomposition is the dynamics resulting of such a quench. Cahn and Hilliard 
have studied the early times of this dynamics by linearizing equation ([1]) around \I/ = (i.e. 
neglecting the non linear term \l/ 3 ). Considering \1/ as a sum of Fourier modes: 

tf(r,t) = ^^ q e iqr+<rf (2) 
q 

where \l/ q is the Fourier coefficient at t = 0, they obtained for the amplification factor er(q) : 

^(q) = -(q 2 + |)q 2 (3) 

This shows immediately that \& = is linearly stable for e > while a band of Fourier 
modes are unstable for negative e, since cr(q) > for < q < \J (s/2). Moreover, the 



most unstable mode is for qc-H = V~ £ /2- This wave number of maximum amplification 
factor will dominate the first stage of the dynamics; this explains in particular why the 
homogeneous domains appear at length scales close to L = \ C -h/2 = n/qc-H, half the 
wave length associated with the instability. For longer times, interfaces separating each 
domain interact through Ostwald ripening, causing < L > to change slowly toward higher 
values Q, Q 

We will now use known results on non-homogeneous solutions of the (G-L) equation to 
study both the saturation of the spinodal decomposition and the coalescence. 



B. Stationary States of the Cahn-Hilliard Dynamics : Symmetric Soliton Lattice 
Solutions 

For e < 0, there exists a whole family of solution of the one dimensional (G-L) equation : 

+ 2^ 3 - V 2 ^ = (4) 

2 

These solutions, the so-called soliton-lattice solutions, are : 



■ x , s .., , A _! Lk 2 + 1 



* fc , e (a;) = fcASn(~, k) with £ = A" = y/ 2— — (5) 

where Sn(x, k) is the Jacobian elliptic function sine-amplitude, or cnoidal mode. This family 
of solutions is parametrized by e and by the Jacobian modulus k £ [0, 1], or "segregation 
parameter" . These solutions describe periodic patterns of period 

f% dt 

X = AK{k)i, where K(k) = / = (6) 

Jo VI - k 2 sin 2 1 



is the complete Jacobian elliptic integral of the first kind. Together with k, it characterizes 
the segregation, defined as the ratio between the size of the homogeneous domains, L = A/2, 
and the width of the interface separating them, 2£. The equation ([6]) and the relation 
£ = A -1 , enable to rewrite this family as : 

T . . AK{k)-k a AKik) 

*k,\(x) = y — Sn(^^x, k). (7) 

As for k = 1, Sn(x, 1) = tanh(x), we recover the usual solution 



^i,eW = ^tanh(^x). (8) 



associated with a single interface (or soliton) of width 2/y |e| that connected the two ho- 
mogenous phases \l/ = Because K{1) diverges, it corresponds to a strong, or macro- 
scopic, segregation. In the opposite limit (weak segregation regime), ^k,e( x ) describes a 
sinusoidal modulation 



*/,._„.,(.(•) = k \ ! — sin(J = /^sin(^x) = kqsin(qx) (9) 



It will correspond to the Fourier mode q — ^ = of the initial white noise, with van- 
ishing amplitude \l/ q = kq. Since experiences, numerical simulations and linear stability 
analysis show that A, the spatial period of the pattern is constant during the whole spinodal 
decomposition process, we choose A to coincide with the most instable wave length obtain 
with the Cahn Hilliard linear approach, A = Xq-h — ^ £Q , where e is the quench tem- 
perature. We then obtain a one parameter family of profiles ty*(x, k) = ^/k,x C - H ( x ) which 
describe very well both the linear growth and the saturation of the spinodal decomposition 



17| | . Thus, the dynamics is now reduced to the time evolution of the single free parameter: 
k(t). Using equations (jSJ) and ©, we find that A, k and e are related to one another through 
the state equation 

e(k) ^ 2(l + e)(^f)\ (1 o) 

This implicit equation tells us that if we impose A = Xc-h, the dynamics can be reduced 
to the evolution of e(t), interpreted as the fictitious temperature or "local temperature" 
associated with the segregation parameter k(t) which characterizes the pattern. This tem- 
perature can be extracted from the profile at any time, using the correspondence between 
e and k of equation ffTUl) . For instance, at t — 0, the amplitude is small and we find that 
k(t = 0) = * q ^~ g -> and thus e*(0) = 8tt 2 /A 2 , different a priori from e (e*(0) = f for 
A = X c ~h) '■ the system is initially out of equilibrium 

Somehow, the dynamics of (C-H) can be projected at first order onto a dynamics along the 
sub-family \I/*(x, k) = ty^\ C H {x), which can be considered as an attractor of the solutions, 
i.e. the density profile of the system will evolve with time, staying always close to a function 
k). Thus, using equation (j3J), the (C-H) dynamics reduces to equation : 

The non linearity has now disappeared into ty*(x,k). Using a solubility condition, it is 
possible to transform this equation into a differential equation for k(t) and to compute the 



full non linear part of this dynamics (i.e. the saturation of the spinodal decomposition), 
which leads the system in a well defined stationary state 



C. Saturations of the Spinodal Decomposition Dynamics 



The spinodal decomposition dynamics will saturate when the fictitious temperature e(t) 
will reach the real thermodynamic one, i.e. the quench temperature 6q; that is, using 
equation of state (ITU]) for A = Xc-h, when k = k$ =0.687 solution of the implicit equation : 



2{l + kf)K{kl 



e \ 2 



C-H 



16 



7T 



(12) 



Using linear stability analysis, Langer has shown that the stationary profile thus obtained, 



^k^x C - H ( x )y i s destroyed by stochastic thermal fluctuations. He has identified 
the most instable mode as an "antiferro" mode, leading to a period doubling. The result of 
this destabilization is another periodic profile of alternate interfaces of period A = 2Xc-h, 
where the length of the domains is : L = A/2 = Xc-h- This means that the new stationary 
profile is given by ^k?,2X -H where k{ =0.985 is solution of the implicit relation 



2(1 + kf)K(kt 



S \2 ~ £ o(2Xc-h) 



16 



4tt 2 



(13) 



Again, this new stationary profile turns out to be linearly instable with respect to an "an- 
tiferro" perturbation of period 4Xc-h- 

Thus these families of profiles and instabilities enable to describe the one dimensional 
coarsening as a cascade of doubling process. Each of these successive intermediate profiles 



can be described by an element of the above family of soliton lattice $ 



kt ,2 n xA 



C-H 



X 



The 



sequence of associated segregation parameters {k s n }, are determined by the implicit relations 



„s\2 



2(1 + k s n 2 )K(k s n 



We have found numerically for the first of them 



e (2 n X 



C-H) 



16 



7T Z . 



(14) 



18| 



fcg = 0.6869795924 



fcf = 0.9851675587 



M= 0.99997210165 



k s 3 = 0.9999999999027 



(15) 
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We see that {k%} converges toward = 1 (single interface case), while the amplitude of the 



modulation converges toward v ^ £ °1 . For large n, we can conclude from the implicit relation 



(HM that the ratio of the domain size to the interface width characterized by K(k^) behaves 
as 7r2 n_1 and {k^) 2 = 1 — 16 exp(— 7r2 n ). Each of the stationary profiles 

e \k s n K(k s n ) y^\K{k s n ) VW\K a < VW\% 



2"tt 2«tt n ' y/2(T+W) \/2(T+k. 

(16)' 

is identically destroyed by the Langer "antiferro" instability . 

III. ID OSTWALD RIPENING 

A. Non-symmetric soliton lattice profile as an ansatz for the ID coarsening process 

In order to describe one step of the coalescence process, i.e. the dynamics that starts 
from ty„.(x) and ends with the profile ^ n+ i(x) , we will use another family of equilibrium 



profiles [19|, of period 2A, solutions of (G-L) equation, which write: 



7, ,s K{k) ~. K (k 1 - k' - 1 + k' )VT^Sn (2x, k) 

ib(x, k) = — + ip (x, k) = — , ' i — — (17) 

where k 2 + k' 2 = 1. Note that < ip(x, k) >= whereas the family ip*(x, k) is solution of : 

V 2 ^* = e(k)^- + 2$* 3 + fi(k) (18) 

with fi(a,k) = (k 2 - 1) (^) 3 and e(k) = (k 2 - 5) (^)' (we have dk) > e ). 

Using Gauss' transformation (or descending Landen transformation [20] ) , we can relate 
the soliton lattice of spatial period 2A (and of modulus k) to the soliton lattice of period A 
(and of modulus /z = j^r) as follows 

${x-^k)+${x+^k) = ^^Sn({4x + X)^,fx) (19) 
$(x - j,k) + $(x + i k) = ^Sn(Ax^ } k) (20) 

Note that K{k) = (1 + fi)K(fi). As if k = k s n+1 , /i = k s n) we see that both the initial 
state \l/„(x) = ^k^,2 n x C - H ( x ) an d the final state \l/„ + i(x) = \I/ fc s 1i 2«+ 1 a C -h( x ) °f a s ^ e P °f 
the coalescence process can be describe, modulo a phase shift, by the same function : 

${x, k, 0) = f (x - (1 - 0/2)^, k) + $(x + (1 - 0/2)^, k) (21) 
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with k = k s n+l and A = 2 n Xc-H (see also Figure 1). Therefore we can describe the coalescence 
by a transformation at constant segregation parameter k, while the degree of freedom 0, 
associated with the relative phase between the two profiles, evolves in time from to 1 (or 



— 1) according to the C-H dynamics 22J 



B. Energy Landscape 

In order to prove the usefulness of this ansatz, we have plot the energy averaged over a 
period, J 7 gl(4 > ) — f Fgl(§(%, k, 4>))dx, as a function of the parameter 0, keeping k constant. 
We see in Figure 2 that the value = corresponds to a local maximum of energy, while 
0=1 (or —1) is a minimum. Note that there is no energy barrier in this particular energy 
landscape, in agreement with linear stability analysis. 

Using the properties of the ansatz, we can extract 0(t) from a numerical simulation. 
Indeed 

$(- k, 0) = $(<f>± k) + -0(A - 0^, k) = 2#^, k) (22) 

K 0) = m\ - A, k) + $(-<f>±, k) = 2$(-<j>± k) 

For = 0, $(|) = $(-f ) = ^(k' - 1) while $(§) = -$(-f ) = for = 1. A k = k s n+l 
is known, from A(t) G [0, 1] defined as follows 

2A(t) _ 1 HhKm - Ak'Snj^k) 

$(-|,fc,0(t)) ^r^ + 2k>Sn(0^,k)-kVTTWSn 2 (<P^,ky 

one can extract <b : 



fc / (A-l) + V / A 2 + fc' 2 -2Afc' 2 

m= mh ^i- t %vi-/ (24) 

One sees in Figure3 that the dynamics is dominated by the time required to leave a 
stationary state. Moreover, the dynamics ends when reaches 1 : we then have $ = ^ n +± 

. ...... - V72,Tf . _ n 

2 n+1 



for which f # n+1 + 2^3 _ v 2 ^ n+1 = 0. 



C. Linear stability analysis 



If we look at the time evolution of the profile $(x, 0), starting from the region = 0, 

| $(x,££ +1 ,0) 

9 



we can transform the (C-H) equation into a phase field equation, replacing 4r $(x, k. 



by ^$(a;, 0(t)) x ^ (with k fixed). The dynamics will then be similar to spinodal 
decomposition (eq. (TTTT) ) . for growing initially as exp(t/r n ) and saturating later at = 1. 

In Figure 4 is plotted ff(x, k^ +1 ) = j^l(x, k^ +1 ) for = 0, which corresponds to the 
most unstable mode founded in Langer's linear stability analysis and is characterized by the 
alternated growth and decrease of domains ("antiferro" mode) : 



* L (x, k) = $(x - (1 - 0/2)^, k) - $'(x + (1 - 0/2)^, k) 



where 



verifies 



— -/2K(k)\ 2 Cn(2x^,k)Dn(2x 



ip'(x, k) 



K(k) 



^ A ' (l- VT^k f Sn(2x 



k) 



K(k) 



,k) 



ip'(k, x) 



+ 6^* 2 i;'(k,x). 



(25) 



(26) 



(27) 



In order to describe the evolution of the phase 0(t) we linearize the Cahn-Hilliard 
dynamics around ^ n (x) = $(x, k s n+1 ,4> = 0) inserting ^ (x, t) = ^ n (x) + §f0(i) = 
*& n (x) + j^l(x, k^ +1 )(p(t) into (pQ). We then have the following dynamics 



d0 
dt 



d 2 



dx 2 V 2 



(28) 



where £ is the Lame operator, 
eigenfunction of period 2X C -h 



ven if this operator doesn't have simple (algebraic) exact 
2l\ , ty L (x,k,<fi), for = and = k s n+1) happens never- 



theless to be a good approximation for the eigenfunction of lowest eigenvalue 22|. Due 
to the concavity of J i 'gl(4 ) ) around = 0, (see Figure 2), this eigenvalue will be negative, 
triggering a linear destabilization of the pattern ty n and an exponential amplification of the 
perturbation, i.e. an exponential growth of the translation with time. 
Indeed, using equation ( TlBl we get 

e(k). 



V 2 ^ L 



So 



e -e(k) 



*l + 6(^;v+-V'-Y-)- 



(29) 



,*2„;./ 



(30) 



+ 6 



,2if(fc) 21T(fc) 
■ A j +A A 
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It turns out that for 6 = 0, 



together with 



So finally 



c 2 fc - rS- - ( 4 - fc ' 2 - 4 *0 (^) 2 (Sh - 9- 



^ L + 6v|>* 2 v|/ L - V 2 ^ ~ -6k I — ^ ) (31) 



2 



As the period associated with is 2 A, equation (fit)]) gives ^1 = — 2(i+P) anc ^ ^hus 
equation (128!) can now be written in a simpler form, similar to equation (TlTI) : 

*"dH°TT^ W a^ < 32) 

We recover that $l is an eigenstate of Lame equation with a negative eigenvalue ■ 
This eigenvalue goes to zero as the coalescence progresses towards higher segregation. 

The characteristic time for one step of period doubling is thus r n ~ ^ 2 ~y 2 (the extra eq/2 
comes from ^^l)- As for small k', k' ~ 4exp(— K(k)), we have the following relationship 
between r" 1 and the period of the profile : 

1 - jeft = 3e 2 eM-K(k)) = 3e 2 exp(-2^ 1 7 r) = 3e 2 exp(--f^). (33) 

As lnr n goes like a power of n, we can conclude that in D = 1, the size of the domain will 
evolve for long t as ^Xc-h hi(3£gt). 



IV. APPLICATION : MODULATED PHASE SYSTEMS 

We can use the preceding ansatz to work out the period of modulated phase systems for 
which there is a competition between two types of interactions: a short-range interaction 
which tends to make the system more homogeneous together with a long-range one, or a 
non-local one, which prefers proliferation of domain walls. This competition results in a 
microphase separation with a preferred length scale [l|. These systems can be study using 
a modified Landau- Ginzburg approach, derived from Cahn-Hilliard equation and often use 



for numerical simulations [231 ]: 



0* = (V 2 ^G^) ) _ ^ = V 2 ( £0^ + 2v]> 3 _ V 2^ } _ ^ (34) 
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— /3 2 ty models in the Cahn-Hilliard equation the long-range interactions, which prevents 
the formation of macroscopic domains and favors the modulation. We could have chosen 
other ways of representing this long-range interaction, but the inclusion of such a term, 
following Oono, enables to describe the behavior of modulated systems at T much lower 
than T c , as we will show below. If we suppose for example that in a 3D problem, the long 
range interaction decreases like the full free energy density writes 

= F GL + F mt = -(W(r)) 2 + U>\r) + -^ 4 (r) + / ^(r')g(r', r)*(r)dr' , (35) 

where g(r', r) = y_ r \ -The long range interaction g(r', r) corresponds to a repulsive inter- 
action when \Er(r') and ^(r) are of the same sign : thus it favor the formation of interphases. 
If we want to study the dynamic of this phase separation, we use the Cahn-Hilliard equation : 



^ = V ' (^) = V ' (f * + 2 *° ~ V ^ + J *V)9V> ■ ( 36 ) 



If one recalls that pjh is the Green's function associated with the Laplacian operator V 2 . 
in 3D, the preceding equation then transforms into 

J y(r')g(r',r)dr'^J = J ^(r')V^(r', r)dr = -(3 2 J V(r')6(r', r)dr = -(3 2 ^(r). 

(37) 

which leads to equation (|34p . The family ffTBT) is not anymore an exact stationary solution of 
this dynamics. Nevertheless, ^ n (x, k n , 2 n \ C -H) is an approximate solution and thus can be 
use as a tool for the calculation using a solubility condition (i.e. if we project the dynamics 
on x G Ker(-^C)). There exists a new sequence {k°} which satisfies (e*(k°) — er ) < X I 
d x 2^/ k o > +/? 2 < x | ${,0 >= and which will characterize the new family of approximate 
stationary configuration. When looking at the linear stability analysis of these solutions, 
equation (CQ) now writes 

The profile will become stable when the eigenvalue of this modified Lame's equation become 
negative, i.e. for r" 1 < 1 . This will take place for 3£q exp(— 2 \c H ) — ft 2 - The dynamics of 
period doubling will thus end with a thermodynamically stable ID microseparated phase of 
period A = -^5— ln(y.) 
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V. DISCUSSION ON THE HYPOTHESIS AND CONCLUSION 



We have shown that the choice of two ansatzs within the soliton-lattice family allows a 
reliable description of the one dimensional dynamics of both spinodal decomposition and 
Ostwald ripening. Contrary to our ansatz relies on the hypothesis that during the first 
stage of the dynamics, the periodicity of the order parameter remains constant, while during 
each steps of the coarsening process, it is the parameter k which remains constant. In a 
sense, there are adiabatic ansatzs : the generation of higher harmonics is governed solely by 
the time evolution of k(t) during the non linear saturation of the spinodal decomposition 
and later on solely by <fi(t) during the Ostwald ripening. The validity of the assumptions 



have been investigated in details and checked numerically 

Our analytic method rely on the assumption that at each step of the dynamics, the system 
can be characterized by a specific spatial period : we need therefore to discuss how this 
approach is relevant to the general case where noise is present. We have noted numerically 
that, for the spinodal decomposition, the average size of the modulation is \c-h, with a 
deviation of less than one percent from the value predicted by the linear theory. It does 
not mean that, in a real system, each domain has a length scale of L = Ac_#/2, but that 
the distribution of the domains' length will be centered around L. The coalescence events 
can be neglected during the initial growth of the amplitude of the modulation : as all 
the eigenvalues of the Lame operator (128]) are then positive, the dynamics remains within 
the ansatz subfamily ty*(x,k). Only after this initial growth has saturated, turns negative 
the lowest eigenvalue ; the coalescence process then starts, and dominates the forthcoming 
dynamics. 

We have computed in this article the characteristic time associated with one step of ideal 
coalescence. By ideal coalescence, we mean a process which breaks as few symmetries as 
possible. In a real system, because of initial fluctuations in the periodicity of the pattern 
selected just after the quench, this instability will concern only region of finite size, where 
it choose a certain sublattice, or a range for (for example, varies from to 1), while it is 
the opposite choice in the neighboring region ( varies from to -1). The global symmetry 
is thus recover on the overall (as in an antiferromagnet). During each step of the process, 
the width of the domains will locally double; but for the system as a whole, due to non 
synchronization between regions, the average length scale will vary continuously. 
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We have also shown that in a modified version of the Cahn-Hilliard dynamics, which take 
into account long range interactions, the computation of this lowest eigenvalue enables to 
compute that the spatial period of the thermodynamically stable modulated phase will be 



-eo 



The use of the solubility technics combined with the choice of an adiabatic ansatz might be 
generalized to the study of other non linear dynamics. For instance, spinodal decomposition 
in superfluid Helium or Bose condensate has been argued to be described by a cubic-quintic 
non linear equation 24j ; in this particular case, one needs to retrieve a relevant soliton-like 
family of solution along which to compute the adiabatic dynamics. The same difficulties 
would arise as well when the method will be adapted to higher space dimensions as there is 
no equivalent to soliton lattice in 2D; the ansatzs ( 1T61) and ( 121"T) could however be used in a 
numerical simulation as tools for following the dynamics. 

It might also be possible to extend this approach to non symmetric profile like the ones 



introduced in different contexts 19 
periodic [2a] or spatially periodic 



251 ] . Or for the cases of special quenches which are time 



27]. 
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VI. FIGURES 



Figure 1 

Construction of the two first steady solutions of the (C-H) dynamics, using a superposition of 
the non-symmetric profile i[)*[k,x), itself stationary solution of the (C-H) equation. By changing 
the phase shift between the two profiles entering into the linear combination, one obtains two 
different symmetric profiles, of periods A and segregation parameter k^ = 0.687 (equation (|19p ) or 
of period 2A and segregation parameter A;f = 0.985 (equation (|2U|) ). 

Figure 2 

Profile of the free energy landscape during a coarsening process, F(0). It starts at = for a 
configuration characterized by the segregation ratio kf = 0.687 for which the energy per unit length 
is F(0)~ —0.135; one sees that in this region, the free energy is a concave function of and thus, 
the associated pattern is linearly instable. The elementary step of the coarsening process ends for 
= 1 associated with a pattern characterized by the segregation ratio k^ = 0.985 for which the 
energy per unit length is F(0)~ —0.45. In the region 0=1, the free energy is a convex function 
of 0. Note that there is no energy barrier. 

Figure 3 

Parameter as a function of time. is extract from a numerical integration of the Cahn-Hilliard 
dynamics using relation (124ft to relate at a given time the profile with the phase 0. It starts at 
= with an exponential growth and saturates at = 1. 

Figure 4 

Langer's most instable perturbation mode of destabilization of the soliton lattice is identified 
with = $(x, k, 0) at = 0. It is composed of two antisymmetric patterns, plotted in dotted 
(plain) line, evolving toward right (left) at velocity +^ ( _ ^)i causing an "antiferro" instability 
leading to a period doubling of the pattern. They are the spatial derivative of the initial non 
symmetric profile ip*(x) which has been used to construct our ansatz in Figure 1. 
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1.2 
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